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Abstract: We calculate the bulk viscosity, drag force and jet quenching parameter 
in Improved Holographic QCD. We find that the bulk viscosity rises near the phase 
transition but does not exceed the shear viscosity. The drag force shows the effects 
of asymptotic freedom both as a function of velocity and temperature. It indicates 
diffusion times of heavy quarks in rough agreement with data. The jet quenching 
parameter values computed via the light-like Wilson loop are in the lower range 
suggested by data. 
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1. Introduction 

A novel window in the physics of the strong interactions has been provided recently 
by the experimental efforts at RHIC, . The consensus on the existing data is 
that shortly after the collision, a ball of quark-gluon plasma (QGP) forms that is 
at thermal equilibrium, and subsequently expands until its temperature falls below 
the QCD transition (or crossover) where it finally hadronizes. Relativistic hydrody- 
namics describes very well the QGP |Q, ^, with a shear- viscosity to entropy density 
ratio close to the universal value suggested by the holographic formulation of = 4 
SYM, i. 

The QGP is at strong coupling, and it necessitates a treatment beyond perturba- 
tive QCD approaches, 0. There are several observables that seem to be important 
in understanding measured features of the collisions. They translate into transport 
properties of the strongly coupled plasma, and reliable methods for the calculation 
are in need. 

A fist class of transport coefficients are viscosity coefficients.^ A general ffuid is 
characterized by two viscosity coefficients, the shear rj and the bulk viscosity (. The 
shear viscosity in strongly coupled theories described by gravity duals was shown to 
be universal, Q. In particular, the ratio rj/s, with s the entropy density, is equal to 

This is correlated to the universality of low-energy scattering of gravitons from 
black-holes. It is also known that deviations from this value can only be generated 
by higher curvature terms that contain the Riemann tensor (as opposed to the Ricci 
tensor of the scalar curvature). In QCD, as the theory is strongly coupled in the 
temperature range Tc < T < 3Tc, we would expect that rj/s Recent lattice 

calculations, agree with this expectations although potential systematic errors in 
lattice calculations of transport coefficients can be large. 

Conformal invariance forces the bulk viscosity to vanish. Therefore the A^ = 4 
SYM plasma, being a conformal ffuid, has vanishing bulk viscosity. QCD on the other 
hand is not a conformal theory. The classical theory is however conformally invariant 
and asymptotic freedom implies that conformal invariance is a good approximation in 
the UV. This would suggest that the bulk viscosity is negligible at large temperatures. 
However it is not expected to be so in the IR: as mentioned earlier lattice data 
indicate that in the relevant RHIC range 1 < ^ < 3 the QGP seems not to be a 

-t c 

^ These are the leading transport coefficients in the derivative expansion. There are subleading 
coefficients that have been calculated recently for A/" = 4 SYM, ||^. However, at the present level 
of accuracy, they cannot affect substantially the comparison to experimental data, H . 
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fully conformal fluid. Therefore the bulk viscosity may play a role near the phase 
transition. 

So far there have been two approaches that have calculated the bulk viscosity in 



YM/QCD, ^, 10, 11] and have both indicated that the bulk viscosity rises near 



the phase transition as naive expectation would suggest. The flrst used the method 
of sum rules in conjunction with input from Lattice thermodynamics, P, P|, pUf . 
It suggested a dramatic rise of the bulk viscosity near Tc although the absolute 
normalization of the result is uncertain. The reason is that this method relies on 
an ansatz for the density associated with stress-tensor two point functions that are 
otherwise unknown. 



The second method |Tl[] relies on a direct computation of the density at low 
frequency of the appropriate stress-tensor two-point function. As this computation 
is necessarily Euclidean, an analytic continuation is necessary. The values at a flnite 
number of discrete Matsubara frequencies are not enough to analytically continue. 
An ansatz for the continuous density is also used here, which presents again a po- 
tentially large systematic uncertainty. 

We will see in the present work that our flndings support a rise of the bulk 
viscosity near but the values are much smaller than previously expected. Studies 
of how this affects hydrodynamics at RHIC, |jl3| suggest that this implies a small 
fall in radial and elliptic flow. 

Another class of interesting experimental observables is associated with quarks, 
and comes under the label of "jet quenching" . Central to this is the expectation that 
an energetic quark will loose energy very fast in the quark-gluon plasma because 
of strong coupling. This has as a side effect that back-to back jets are suppressed. 
Moreover if a pair of energetic quarks is generated near the plasma boundary then 
one will exit fast the plasma and register as an energetic jet, while the other will 
thermalize and its identity will disappear. This has been clearly observed at RHIC 
and used to study the energy loss of quarks in the quark-gluon plasma. 

Heavy quarks are of extra importance, as their mass masks some low-energy 
strong interaction effects, and can be therefore cleaner probes of plasma energy loss. 
There are important electron observables at RHIC, that can probe heavy-quark 
energy loss in the strongly coupled quark-gluon plasma. Such observables are also 
expected to play an important role in LHC |T5[ . 

A perturbative QCD approach to calculate the energy loss of a heavy quark in 
the plasma has been pursued by calculating radiative energy loss, ||l^. However its 
application to the RHIC plasma has recently raised problems, based on comparison 
with data. A phenomenological coefficient used in such cases is known as the jet 
quenching coefficient g, and is defined as the rate of change of the average value of 
transverse momentum square of a probe. Current fits, |jl^, |1^ , indicate that a value 
of order 10 GeV'^/fm or more is needed to describe the data while perturbative 
approaches are trustworthy at much lower values. 
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Several attempts were made to compute quark energy loss in the holographic 
context, relevant for A/" = 4 SYM^. In some of them [|19], the jet-quenching 



coefficient q was calculated via its relationship to a light-like Wilson loop. Holography 
was then used to calculate the appropriate Wilson loop. The q obtained scales as 
\/A and as the third power of the temperature, 

^ [i] .ATT 



^conformal — „ 1-5-1 v2A Tt'^T (1.1) 

uJ 

A different approach chooses to compute the drag force acting a string whose UV 
end-point (representing an infinitely heavy quark) is forced to move with constant 
velocity f , 1^1], 0, in the context of A/" = 4 SYM plasma. The result for the drag 
force is 

-^conformal = TT"^ = (1-2) 



2 



and is calculated by first studying the equilibrium configuration of the appropriate 
string world-sheet string and then calculating the momentum fiowing down the string. 
This can be the starting point of a Langevin evolution system, as the process of 



energy loss has a stochastic character, as was first pointed out in [H and more 
recently pursued in p5l-pT|. 



Such a system involves a classical force, that in this case is the drag force, and 
a stochastic noise that is taken to be Gaussian and which is characterized by a 
diffusion coefficient. There are two ingredients here that are novel. The first is that 
the Langevin evolution must be relativistic, as the quarks can be very energetic. Such 
relativistic systems have been described in the mathematical physics literature, 
and have been used in phenomenological analyses of heavy- ion data, ||l3 . They are 



known however to have peculiar behavior, since demanding an equilibrium relativistic 
Boltzmann distribution, provides an Einstein relation that is pathological at large 
temperatures. Second, the transverse and longitudinal diffusion coefficients are not 
the same, [Q. A first derivation of such Langevin dynamics from holography was 



given in [^. This has been extended in in where the thermal- like noise was 
associated and interpreted in terms of the world-sheet horizon that develops on the 
probe string. 

Most of the transport properties mentioned above have been successfully com- 
puted in A/" = 4 SYM and a lot of debate is still waged as to how they can be applied 



to QCD in the appropriate temperature range, [Q,|3^,|^. A holographic descrip- 



tion of QCD has been elusive, and the best we have so far have been simple bottom 
up models. 

In the simplest bottom- up holographic model known as AdS/QCD the bulk 
viscosity is zero as conformal invariance is essentially not broken (the stress tensor 



Most are reviewed in |18 
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is traceless), and the drag force and jet quenching essentially retain their conformal 
values. 



In the soft-wall model no reliable calculation can be done for glue correlators 
and therefore transport coefficients are ill-defined, as bulk equations of motion are not 
respected. Similar remarks hold for other phenomenologically interesting observables 
as the drag force and the jet quenching parameter. 

A hybrid approach has been advocated in ^ 5^ combining features of 



bottom-up and top-down (string theory) models. Such an approach is essentially a 
five-dimensional dilaton-gravity system with a non-trivial dilaton potential. Flavor 
can be eventually added in the form of Nf space-time filling D4 — D4 brane pairs, 
supporting U^Nf)^ x U{Nf)ji gauge fields and a bi-fundamental scalar [pH '^. 



The UV asymptotics of the potential are fixed by QCD perturbation theory, 
while the IR asymptotics of the potential can be fixed by confinement and linear 
glueball asymptotics. 

An analysis of the finite temperature behavior M^, HH] has shown that the phase 



structure is exactly what one would expect from large- A^^c YM^. Einstein-dilaton 
gravity with a strictly monotonic dilaton potential that grows sufficiently fast, generi- 
cally shares the same phase structure and thermodynamics of finite-temperature pure 
Yang-Mills theory at large N^. There is a deconfinement phase transition (dual to 
a Hawking-Page phase transition between a black-hole and thermal gas background 
on the gravity side), which is generically first order. The latent heat scales as A"^. 
In the deconfined gluon-plasma phase, the free energy slowly approaches that of a 
free gluon gas at high temperature, and the speed of sound starts from a small value 
at Tc and approaches the conformal value cj. = 1/3 as the temperature increases. 
The deviation from conformal invariance is strongest at T^, and is signaled by the 
presence of a non-trivial gluon condensate, which on the gravity side emerges as a 
deviation of the scalar solution that behaves asymptotically as close to the UV 
boundary. In the CP-violating sector, the topological vacuum density tr FF has 



zero expectation value in the deconfined phase, in agreement with lattice results 
and large- A^c expectations. 



The analysis performed in was completely general and did not rely on any 



specific form of the dilaton potential V^(A). A potential with two parameters, was 
subsequently chosen to describe YM data, The (dimensionless) free energy, 

entropy density, latent heat and speed of sound, obtained on the gravity side by nu- 
merical integration of the 5D field equations, were compared with the corresponding 
quantities, calculated on the lattice for pure Yang-Mills at finite-T, resulting in ex- 
cellent agreement, for the temperature range that is accessible by lattice techniques. 
The same model also shows a good agreement with the lattice calculation of glue- 



^D4 — D4 brane pairs for flavor where first suggested in and the finite temperature sohitions 
studied in [||. 

^Similar resuhs, but with somewhat different potentials were also obtained in Efl, B7|. 
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ball mass ratios at zero temperature. Moreover the value of the deconfining critical 
temperature (in units of the lowest glueball mass) was also in good agreement with 
the lattice results. 

In short, the model, named Improved Holographic QCD, (or IHQCD for short), 
gives a good phenomenological (holographic) description of most static properties 
(spectrum and equilibrium thermodynamics) of large- iY^ pure Yang-Mills, as com- 
puted on the lattice, for energies up to several times T^. Therefore it constitutes 
a good starting point for the computation of dynamical observables in a realistic 
holographic dual to QCD (as opposed to e.g. Af = 4 SYM), such as transport coef- 
ficients and other hydrodynamic properties that are not easily accessible by lattice 
techniques, at energies and temperatures relevant for relativistic heavy-ion collision 
experiments. 

The purpose of the present paper is to compute transport properties (the bulk 
viscosity) and energy loss coefficients (the jet quenching parameter and the drag 
force) in the specific Improved Holographic QCD model described in . 



The shear viscosity of IHQCD is the same as that of A/" = 4 SYM, as the model 
is a two derivative model. Although this is not a good approximation in the UV of 
QCD, it is expected to be a good approximation in the energy range < T < 5Tc. 
We find that the bulk viscosity rises near the phase transition but ultimately stays 
slightly below the shear viscosity. We also give a general holographic argument that 
any (large-N) gauge theory that confines color at zero temperature should have an 
increase in the bulk viscosity-to-entropy density ratio close to Tg. 

The drag force on heavy quarks, and the associated diffusion times, are calcu- 
lated and found to be momentum depended as anticipated from asymptotic freedom. 
Numerical values of diffusion times are in the region dictated by phenomenological 
analysis of heavy-ion data. We calculated the medium-induced corrections to the 
quark mass (needed for the diffusion time calculation), and we find they result in 
a mildly decreasing effective quark mass as a function of temperature. This is con- 
sistent with lattice results. Finally the jet-quenching parameter is calculated and 
found to be comparable at Tc to the one obtained by extrapolation from A/" = 4 
SYM. Its temperature dependence is however different and again reflects the effects 
of asymptotic freedom. 

There are several sources of error and systematic uncertainties in the results 
above. We analyze them in the appropriate sections and make a long commentary 
on their importance in the last section. 

The structure of the paper is as follows. In Section 2, we review the holographic 
construction that shall be used to compute certain observables of QGP. In particu- 
lar we review the asymptotic behaviors of the backgrounds and discuss the various 
parameters in the model and how they are flxed. Section 3 is devoted to the com- 
putation of the bulk viscosity. We describe the general holographic computation of 
the quantity from the graviton fluctuation equations on the dual background and 
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compute its temperature dependence numerically. We also make a proposal for a 
holographic explanation of the rise in the bulk viscosity near the phase transition. 
In section 4, we compute the drag force on a heavy quark moving in QGP in our set- 
up. We obtain general analytic formulas in the relativistic and the non-relativist ic 
limits of the drag force as a function of temperature. We compare our findings with 
the A/" = 4 SYM result. In this section we also compute the diffusion times for the 
heavy quarks in the QGP. In particular, we present numerical results for the charm 
and the bottom quarks. Furthermore, we compute the thermal corrections to the 
quark masses in our set-up and finally discuss in what temperature ranges should 
our results be trusted. In section 5, we compute the jet-quenching parameter in our 
set-up. Our findings are compared with the conformal (A/" = 4 SYM) case. Finally, 
the section 6 contains a discussion and outlook. The various appendices detail our 
computations. 



2. Review of IHQCD backgrounds 



The holographic duals of large Nc Yang Mills theory proposed in ||38|, |3^ are based 
on five- dimensional Einstein-dilaton gravity with a dilaton potential. The basic fields 
for the pure gauge sector are the 5D metric (^^i, (dual to the 4D stress tensor) and a 
scalar field $ (dual to TrF"^). The Einstein frame action for these fields is: 

^5 = -M^N^ [ d^'x^ R - \{d^f + V{^) + 2MlNl ! d^xVh K. (2.1) 

J I J JdM 

Here, Mp is the five-dimensional Planck scale and is the number of colors. The 
last term is the Gibbons-Hawking term, with K being the extrinsic curvature on the 
boundary. The effective five-dimensional Newton constant is = l/llGirMpN^), 
and it is small in the large- iVc limit. 

The scalar potential V{^) is what determines the dynamics. Its form is in part 
motivated from non-critical string theory, and in part chosen following guidelines 
from phenomenology. We will often write ^ as a function of A = e*. 

Asymptotic freedom in the UV requires V{X) have a regular expansion for small 

' " 12 

V(X) = -^{l + V„\ + ViX' + ...), Vo>0, A^O. (2.2) 

This ensures that in any solution of Einstein's equations the metric has an asymp- 
totically AdS^ UV region, with AdS length i, in which the field A vanishes logarith- 
mically. We have the perturbative identification, valid for small A: 

k\ = N^gyM, = (2-3) 

o Po 

where Pq = (22/3) (47r)^^ is the first coefficient of the perturbative beta-function for 
the 't Hooft coupling NcQym pure Yang-Mills theory. 



-7- 



For large A, confinement and a linear glueball spectrum require: 

V{X)^V^X'/^{\ogXy/' A^oo, (2.4) 

where Voo is a positive constant^. 

For a generic potential that obeys the asymptotics (|2.2| ) and (|2.4|) , the model 
exhibits the following features: 

• Zero temperature. At zero temperature, the gravity solution is dual to a 
confining 4D theory. One has color confinement (i.e. a Wilson Loop area law) 
and a discrete glueball spectrum with linear behavior, ~ n. The solution 
of Einstein's equations for the metric and dilaton has the form: 

ds^ = e^^-^") {dr^ + r^^^dx^'dx") , A = Ao(r), < r < oo, (2.5) 

with small-r log-corrected AdS asymptotics, 

A,Xr) ~ - log rll + O (^) + . . . , A„(r) ~ (2.6) 



and large-r behavior: 

^3 



~ -Cr , Ao(r) ~ exp 



2 



r ^ oo. (2.7) 



In equation (|2.6|) , A is an integration constant that sets the length scale of 
nonperturbative physics; the constant C in equation ( p7|) is determined in 
terms of A. 

Finite temperature. At finite temperature, one finds a first order phase 
transition between a low-temperature confined phase, described by the solution 



(|2.5|) , and a high-temperature deconfined phase, described holographically by 
D black hole 

ds^ = e^^M 



a 5D black hole solution: 

dr"^ 



fir) 



- f{r)dr + dx^dx"^ , A = A(r), < r < r^. 

(2.8) 

These solutions are characterized by the presence of a horizon where f{rh) = 
0, and have a temperature T and an entropy density s: 

T = -l^^ s = 47r(MX)e'^^''^. (2.9) 

In the UV (r — )• 0), and for any Vh , the black holes are asymptotically AdS^ 
and reduce to the zero-temperature metric A{r) ^ Ao{r), f{r) ~ 1. 



^Other types of large-A asymptotics also lead to color confinement, with different features of the 
^lueball spectrum. These solutions were analyzed in complete generality in [p9[. 
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In all types of solutions, (|2.5|) and ( p.8|) , the dilaton $(r) is a monotonically 
increasing function of r. One can therefore use $ itself as the radial coordinate in 

ds^ = e^M^\-fde + dxmdx"^) + e^^W^. (2.10) 



Comparison of ( p^ ) and ( |2.10| ) determines, 

E = A - lo£ 



(2.11) 



This form of the metric will prove useful later. 

Generically, in these types of models there exist two separate black-hole solutions. 



that were referred in as the hig and the small black-holes. In it was proved 
that existence of this second type of black-hole solution (the small BH) is necessary 
and sufficient for a first order confinement-deconfinement phase transition. 

The big BH solution exists for T > T^m for some finite T^m, see figure ^ (b), 
and dominates the entire thermodynamic ensemble for T > where is always 
larger than T^j^. It always dominates in the thermodynamic ensemble over the small 
BH. It corresponds to the range < < Xmin in the horizon value of the dilaton, 
for some finite Xmin, see figure ^ (b). This solution is proposed as the holographic 
dual of the Yang-Mills gluon plasma. 

The small BH solution also exists for T > Tmin and it corresponds to the range 
Amin < Aft < oo, see figure ^ (b). As it is never dominant in the thermodynamic 
ensemble, it bears no direct significance for an holographic investigation of the quark- 
gluon plasma. 

This situation is depicted in Figure |l|. 

In summary, there exists three separate solutions to the dilaton-gravity system: 



i. The thermal gas ( p.5|) that exists for all T > 0. It is the dominant solution for 

T < Tr. 



ii. The big BH (|2.8|) that exists for T > T^m and becomes the dominant solution 
for T > T,. 

iii. The small BH that exists for T > Tmin and is always sub-leading in the ther- 
modynamic ensemble. 



The solutions ( |2.5| ) and ( |2.8| ) are written in the Einstein frame. Some of the 
transport properties we compute in this paper however are defined in terms of the 
string frame, since they are related to world-sheet quantities. In the five-dimensional 
non-critical string setup, the string frame and the Einstein frame metrics are related 
by 0: 

dsl = et*d4, (2.12) 
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F 



Nl 




Figure 1: Free energy of black hole solutions as a function of temperature. The (constant) 
free energy of the T = confining vacuum is set to zero. The two branches correspond 
to the big black holes (lower branch) and the small black holes (upper branch). The two 
branches merge at a minimum temperature Tmin > 0, corresponding to the vertical dashed 
line. The free energy of the big black hole branch crosses the x-axis at T = Tc, indicating 
a first order phase transition between the vacuum and big black hole phase. 



so we can define a string frame scale factor (both at zero and finite temperature): 

Mr) = Air) + ^$(r). (2.13) 



As shown in the Einstein frame scale factor is monotonic if the metric is 



asymptotically AdS and the theory satisfies the null energy condition^. On the other 
hand, the string frame scale factor may not be monotonic. In particular, in the 
backgrounds with IR asymptotics (|2.7| ), (which follows if the dilaton potential obeys 
the zero-temperature string frame scale factor behaves as: 



. , . J - logr/£ +00 r 0, i 
log r — !■ +00 r — 7- oo. 



Therefore, the zero-temperature string frame scale scale factor must have a minimum 
at some finite value of the radial coordinate, r = r*, where in string units the metric 
has a minimum size e'^" "^''*-'. This is what causes the holographic Wilson loop to 



exhibit an area law jBOI. The confining string tension ac is related to the fundamental 



string length is and value of the string frame metric at the extremum: 

Notice that it is not guaranteed that the minimum survives in the black hole 
solutions. For sufficiently high temperature, the minimum of the string world-sheet 



^This is always the case for a single scalar field with a canonical kinetic term. 
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should disappear behind the horizon. In fact, this is what happens in the explicit case 
we will consider in this paper: a numerical analysis shows that for all temperatures 
larger than the critical temperature T^, both the string and Einstein frame scale 
factor are monotonically decreasing over the whole range < r < r/j. Thus the 
minimum of the scale factor in both frames occurs at r = Th- 



in 1^ we assumed a specific form of the potential: 

V^(A) = ^ {l + K,A + V^X'/' [log (1 + + VsX')] . (2.16) 

The model specified by the potential ( p. 161) contains a few adjustable parameters, 
namely the coefficients Vi and i entering the potential, and the 5D Planck scale Mp. 
They were fixed in as follows: 

• The coefficients Vq and V2 are chosen to reproduce the perturbative Yang- Mills 
beta-function up to 2-loop order, /3(A) = — /3oA^ — f3iX^ + O(A^). This requires: 

V„.% V,.pi(?l±^)\ (2,7) 



9' - ' ^ \^ 81^^ 

For pure Yang-Mills the beta-function coefficients are: 

The coefficients Vi and V3 were fixed by comparing the latent heat of the phase 
transition, and the pressure of the deconfined phase at a given temperature 
(T = 2Tc), to the corresponding lattice results. A successful matching leads to 
the choice: 

Vi = 14, V3 = 170. (2.19) 

The asymptotic AdS scale i only affects the overall unit of energy, and can be 
set by fixing the value of a single dimensionful quantity in the model (say the 
lowest glueball mass, or the critical temperature). Any physical dimensionless 
quantity is independent of i. Once i is given, the UV solution is asymptotically: 



' log log r 

r ' " \\ogr J ' ^ /^ologAr ' " \^ log^r 



A(r) = log3 + 0(^), A(r) = -^-f— ^ + 0(^^^). (2.20) 



The scale A appearing in the UV asymptotics of A(r) is an integration con- 
stant of the zero-temperature Einstein's equations, and it is related to the UV 
boundary conditions [Aq, Aq) at a small but finite coordinate Vq as: 



A ^ r\PoXo)-'exp \ Ao - -^1 , b=^. (2.21) 
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It may seem from this discussion that there is an extra dimensionless parameter 
in our model, A^, with respect to pure 4D Yang-Mills (where the only parameter 
is the scale A). This is not so: all physical quantities that can be related 
holographically to a Yang-Mills observable have a trivial dependence on M. In 



fact, as shown in [^, changing A while keeping £ fixed, is the same as shifting 
A{r) by a constant, i.e. a fixed rescaling of all energies in the model or a change 
of units. Alternatively, for any given value of there exists a unique solution 
such that the scale A is equal to the physical value in 4D Yang-Mills, and that 
no dimensionless observable depends on this choice. 

• The 5D Planck scale is fixed (in units of €) so that, in the T oo limit, the 
equation of state matches that of a free relativistic gas of photons, 

lim ^ = !^ ^ M3 = il!-. (2.22) 

As shown in with these choices of the parameters the 5D holographic model 
is able to accurately reproduce all known thermodynamic properties of finite tem- 
perature Yang-Mills theory, as they emerge from lattice studies. It also displays a 
glueball spectrum which is in good agreement with lattice results. The value of the 
confinement- deconfinement transition is found to be = 247 MeV, very close to the 
lattice determination of the YM critical temperature. 

In the following sections we discuss the transport coefficients (i.e. the bulk 
viscosity) of the deconfined phase and the energy loss of a heavy quark in this specific 
holographic model. 



3. Bulk viscosity 



The bulk viscosity ( is an important probe of the quark-gluon plasma. Its profile as 
a function of T reveals information regarding the dynamics of the phase transition. 
In particular, both from the low energy theorems and lattice studies ^, |rT| , there 
is evidence that ( increases near Tc. 

For a viscous fiuid the shear t] and bulk ( viscosities are defined via the rate of 
entropy production as 

2 



ds 
di 



1 
T 



diVj + d^v 



(3.1) 



Therefore, in a holographic setup, the bulk viscosity can be defined as the re- 
sponse of the diagonal spatial components of the stress-energy tensor to a small 
fiuctuation of the metric. It can be directly related to the retarded Green's function 
of the stress-energy tensor by Kubo's linear response theory: 



(3.2) 
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where Gr{w,'p) is the Fourier transform of retarded Green's function of the stress- 
energy tensor: 

3 

Gr{w,p) = -i d^xdte''''-'P^--'e{t)J2{[Tnit,x),T,^^^^ (3.3) 

A direct computation of the RHS on the lattice is non-trivial as it requires analytic 
continuation to Lorentzian space-time. In refs. the low energy theorems of 

QCD, as well as (equilibrium) lattice data at finite temperature were used in order to 
evaluate a particular moment of the spectral density of the relevant correlator, using 
a parametrization of the spectral density via two time-dependent constants, one of 
which is the bulk viscosity a relation for their product was obtained as a function of 
temperature. This can be converted to a relation for assuming the other constant 
varies slowly with temperature. 

The conclusion was that C,/s increases near T^. Another conclusion is that the 
fermionic contributions to C, are small compared to the glue contributions. 

The weak point of the approach of is that it requires an ansatz on the 
spectrum of energy fluctuations, and further assumptions on the other parameters, 
which are not derived from first principles. 



A direct lattice study of the bulk viscosity was also made in |TT|]. Here, the result 
is also qualitatively similar ^ However, the systematic errors in this computation 
are large especially near Tc, mostly due to the analytic continuation that one has to 
perform after computing the Euclidean correlator on the lattice. 

The results of references |B[,H] and the assumptions of the lattice calculation 



have been recently challenged in [|12 . 



3.1 Holographic computation and main results 

The holographic approach offers a new way of computing the bulk viscosity. In 
the holographic set-up, C is obtained from (p.2|). Using the standard AdS/CFT 



prescription, the two point-function of the energy-momentum tensor can be read off 
from the asymptotic behavior of the metric perturbations 6g^y. This is similar in 



spirit to the holographic computation of the shear viscosity |5T[], but it is technically 
more involved. For a recent treatment of the fluctuation equation governing the 
scalar mode of a general Einstein-Dilaton system, see |^ . Here, we shall follow the 
method proposed by ||53| . 

As explained in , one only needs to examine the equations of motion in the 
gauge r = (j), where the radial coordinate is equal to the dilaton. In our type of 
metrics, the applicability of this method requires some clarifications, that we provide 
in App. 0. Using SO (3) invariance and the five remaining gauge degrees of freedom 
the metric perturbations can be diagonalized as 

6g = diag{goo, gn, 911,911,955), (3.4) 
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where 



2B 



900 = -e^^/[l+/ioo(0)e-*"*], gu = e^^[l+/in(0)e-n, 9^5 = —[l+h,{<j))e 

(3.5) 

where the functions A and B are defined in ( p.lO| ). Here, the fiuctuations are taken 
to be harmonic functions of t while having an arbitrary dependence on (p. 

The bulk viscosity depends only on the correlator of the diagonal components of 
the metric and so it suffices to look for the asymptotics of hu. Interestingly, in the 
r = (f) gauge this decouples from the other components of the metric and satisfies the 
following equation'^ 

Ki- f- A - 4^' + ^B' - 41 h\,- i-'-^uj' + 7^ - ^1 /^ii = . (3.6) 



One needs to impose two boundary conditions. First, we require that only the 
infalling condition survives at the horizon: 

/^ll ^ C6(0,, -0)"^, (3.7) 

where Cf, is a normalization factor. The second boundary condition is that hu has 
unit normalization on the boundary: 

/ill ^ 1, (3.8) 

Having solved for /;,ii(0), Kubo's formula (|3.2|) and a wise use of the AdS/CFT 



prescription to compute the stress-energy correlation function |Q determines the 
ratio of bulk viscosity as follows. 

The AdS / CFT prescription relates the imaginary part of the retarded Ta Green 



function to the number fiux of the hu gravitons J-" |]53 



/mG^(u;,0) = --^ (3.9) 

where the fiux can be calculated as the Noether current associated to the U{1) 
symmetry hu — )■ e^^hu in the gravitational action for fiuctuations. One finds, 

J" = ^^Jr^[huh\, - huhU (3.10) 

As J-" is independent of the radial variable, one can compute it at any 0, most 
easily near the horizon, where hu takes the form ( |3.7| ). Using also the fact that 
{dA/d(f)){(j)h) = -8V{(f)h)/9V'{(j)h) (see Appendix 0), one finds 



'Difference in the various numerical factors in this equation w.r.t |53 is due to our different 
normalization of the dilaton kinetic term. 
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Then, (|3.2|) and ( p.9|) determine the ratio of bulk viscosity and the entropy density 
as, 

C_ 3 ( V'i<Ph) \\^,2 
s 327r V VM 



Cbp. (3.12) 



In the derivation we use the Bekenstein-Hawking formula for the entropy density, 
s = exp3A{(f)h)/4:G5. 

To find ( we need to find Cf, only in the limit a; — ?■ 0. The computation is 



performed by numerically solving equation (|3.6| ) with the appropriate boundary con- 
ditions. There are two separate methods that one can employ to determine the 
quantity Cb. 



1. One can solve ( p. 61) numerically with a fixed u/T, but small enough so that Cb 
reaches a fixed value. The method is valid also for finite values of u. From a 
practical point of view, it is easier to solve ( p.6|) with the boundary condition 
( |3.7| ) with a unit normalization factor, read off the value on the boundary 

oo) from the solution and finally use the symmetry of ( |3.(j| ) under constant 
scalings of hu to determine \cb\ = l/|/iii(— oo)|. 

2. An alternative method of computation that directly extracts the information 
at = follows from the following trick . Instead of solving ( p. 6]) for small 
but finite u, one can instead solve it for = 0. This is a simpler equation, yet 
complicated enough to still evade analytic solution. Let us call this solution /i^i- 
One numerically solves it by fixing the boundary conditions on the boundary: 
hii{—oo) = 1 and the derivative dhl^/ dcf) {—oo) is chosen such that hu is regular 
at the horizon. Matching this solution to the expansion of { \i.7\j for small co 
than yields \cb\ = h1^{(j)h). 

We used both methods to obtain (/s as a function of T and checked that they 
yield the same result. As explained in [|^, most of the thermodynamic observables 
are easily computed using the method of scalar variables. This method is summarized 
in appendix ^ where we also detail the computation of the bulk viscosity using these 
variables. 

Our results are presented in figure ^ This figure gives a comparison of the 
curve obtained by the holographic calculation sketched above by solving ( p.6| ) and 
the lattice data of |]rT|. We also show t]/s = l/An in this figure for comparison. The 



result is qualitatively similar to the lattice result where (/s increases as T approaches 
Tc, however the rate of increase is slower than the lattice. As a result, we obtain a 



value (/s{Tc) ~ 0.06 that is an order of magnitude smaller than the lattice result [p 
which is 0.8. Note however that the error bars in the lattice evaluation are large near 
Tc and do not include all possible systematic errors from the analytic continuation. 

We should note the fact that the holographic calculation gives a smaller value 
for the bulk viscosity near than the lattice calculation is generic and has been 
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Figure 2: Plot of C/s (continuous line) calculated in Improved Holographic QCD model. 
This is compared with the lattice data of that are shown as boxes. The horizontal 
dashed line is indicating the (universal) value of - for comparison. 



found for other potentials with similar IR asymptotics, |53|. The fact that the value 
of (/s near is correlated with the IR asymptotics of the potential will be shown 
further below. 

Another fact that one observes from figure ^ is that (/s vanishes in the high T 
limit. This reflects the conformal invariance in the UV and can be shown analytically 
as follows. C/'^ is determined by formula ( |3.12| ). In the high T limit, (corresponding 
to A/i — 7- 0, near the boundary), the fluctuation coefficient \cb\ — ?■ 1. This is because 
of the boundary condition /in (A = 0) = 1. We use the relation between T and A/j in 
the high T limit 

A;,^ (/3olog(7rT/A))-\ (3.13) 
Substitution in ( p.l2| ) leads to the result. 



big 



27n log"(7rr/A) 



as T — )■ oo. 



(3.14) 



As s itself diverges as in this limit - it corresponds to an ideal gas - we learn 
that C also diverges as TVlog^(T). Diver gence at high T is expected from the bulk- 
viscosity of an ideal gas. We do not expect however the details of the asymptotic 
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(a) 



(b) 



Figure 3: (a) Numerical evaluation oi Q/r] both on the big-BH branch (the solid curve) 
and on the small BH branch (the dashed curve). Tm denotes Tmin- (b) The two branches 
of black-hole solutions, that correspond to different ranges of X^. The big BH corresponds 
to \h < Xmin and the small BH corresponds to Xh > Xmin- 

result to match with the pQCD result, for the same reasons that the shear-viscosity- 
to-entropy ratio does not, |Q. However, the T-dependence is very similar to the 



3.2 Holographic explanation for the rise of (^/s near Tc and the small black- 
hole branch 

With the same numerical methods, one can also compute the ratio (/s on the small 
black-hole branch. As this solution has a smaller value of the action than the big 
black-hole solution, it is a subleading saddle point in the phase space of the theory, 
hence bears no direct significance for an holographic investigation of the quark- 
gluon plasma. However, as we show below, the existence of this branch provides a 
holographic explanation for the peak in C/s in the quark-gluon plasma, near Tc. 

From the practical point of view, we find the second numerical method above 
(solving the fluctuation equation at w = 0) easier in the range of Xh that corresponds 
to the small black- hole. The result is shown in figure |] (a). The presence of two 
branches for T > Tmm, is made clear in this figure. See also fig H (b) for the respective 
ranges of Xh that correspond to small and big BHs. In fig |] (a), (/s on the big BH 
branch is depicted with a solid curve and the small BH branch is depicted with a 
dashed curve. We observe that (/s keeps increasing on the big-BH branch as T is 
lowered, up to the temperature T^m where the small and big BH branches merge*. 
On the other hand, on the small BH branch (/s keeps increasing as the T is increased, 

^As far as the thermodynamics of the gluon plasma is concerned, the temperatures below (on 
the big BH branch) has little importance, because for T < T^. the plasma is in the confined phase. 



pQCD result, |5l: 




(3.15) 
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Figure 4: An inset from tlie figure |3| around the maximum of C/s. 

up to a certain T^ax that lies between Tmin and T^, see figure ^. From this point 
onwards, (/s decreases with increasing T. 

A simple fact that can be proved analytically is that the derivative of(/s diverges 
at Tmin- This is also clear from figure |. This is shown by inspecting equation ( |3.12| ). 



The T derivative is determined as d/dT = {dT / d\h)d / dX^. Whereas the derivative 
w.r.t \h is everywhere smooth^, the factor dT/dXh diverges at Tmin by definition, see 
figure H (b). 

Therefore, we propose that the presence of a Tmin where the big and the small 
black-holes meet, in other words, presence of a small-black- hole branch is responsible 
for the increase of (/s near Tmin- As in most of the holographic constructions that 
we analyzed, and specifically in the example we present in this paper, Tc and Tmin 
are close to one another, this fact implies a rise in the bulk viscosity near Tc. This 
proposal, combined with the fact that the existence of a small BH branch and color 
confinement in the dual gauge theory at zero T are in one-to-one correspondence [^ , 
suggests that in confining large-N gauge theories, there will he a peak in the ratio (/s 
close to Tc- 

Another fact that can be shown analytically is that (/s asymptotes to a finite 
value as T — )■ oo in the small black- hole branch We find that. 



— )■ — , as T — )• oo. (3.16) 

small 



As the entropy density vanishes in this limit [^, we conclude that C should vanish 
with the same rate. 



^Note that ct is also a function of A/j. As both the fluctuation equation (3.6) and the boundary 
conditions are smooth at Xh = A,„i.„, one concludes that Cb also is smooth at this point. 
^°See equations (B.l), (B.2), and the discussion in App. B. 
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For a general potential with strong coupling asymptotics 



V{X) ~ as X^oo, (3.17) 

taking into account ( |3.12|) , equation ( p.l6|) is modified to 



"■^ '^'^ ~^ '^0- (3.18) 

small "^^^ 



where tq is the position of the singularity in the zero temperature solution. 

For confining theories, the limit Th — )■ corresponds to T — )■ oo on the small BH 
branch. However, one can show that the result ( |3.18D holds quite generally, regardless 
of whether the zero T theory confines or not^^. In particular, for the non-confining 
theories — that is either when Q < 4/3 or when (5 = 4/3 but the sub leading term 
in the potential vanishes at the singularity — there is only the big black-hole branch 
and the limit — J- tq corresponds to the zero T limit of this BH. Thus, we also learn 
that there exist holographic models that correspond to non-confining gauge theories 
whose zero T limit yield a constant (/s. This constant approaches zero as Q — )■ 0, 
i.e. in the limiting AdS case. 

We also see that the asymptotic value of C/s in the small BH branch is close 
to the value of C/s near T^. We shall give an explanation of this fact in the next 
subsection. Using the asymptotic formula (|3.18|) , the fact that Q > | for confinement 



and Q < for the IR singularity to be good and repulsive we may obtain a range 
of values where we expect C/s to vary, namely 



± < £ 

Gtt s 



< i-. (3.19) 

small, asymptotic 



A final observation concerns the coefficient Cb{Xh) in (p. 12 ). This part is the only 



input from the solution of the fluctuation equation, the rest of (|3.12|) is fixed by the 
dilaton potential entirely. We plot the numerical result for Cb in fig ^ as a function 
of the coupling at the horizon Xh- 

First of all. Figure ^ provides a check that, the approximate bound of \cb\ > 1, 



is satisfied in the entire range. One also observes Cb approaches to 1 in the IR and 
UV asymptotics. These facts can be understood analytically: In the UV (near the 
boundary) it is because of the boundary condition Cb = 1. In the IR, it is more 
subtle, and we explain this in appendix 

Finally, we observe that the deviation of Cb from the asymptotic value 1 is max- 
imum around the phase transition point Ac. In fact, we numerically observed that 
the top of the curve in figure ^ coincides with Ac to a very high accuracy. Whether 
this is just a coincidence or not, it needs to be clarified. 
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The arguments in App. B remain valid in the general case. 
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Figure 5: The coefficient \cb\ of equation ( 3.12 ) as a function of Xh- 



3.3 The adiabatic approximation 



Motivated by the Chambhn-ReaU solutions Gubser et al. proposed an 

approximate adiabatic formula for the speed of sound. In the case when V' /V is 



a slowly varying function of 0, ||5g] proposes the following formulae for the entropy 
density and the temperature: 

8 Z"^" V 

hgs = -- #— + ■■■, (3.20) 



3.1 

/IV' 8V\ 

logr=/ <i^(-----)..., (3.21) 



where ellipsis denote contributions slowly varying in 



12 



It is very useful to reformulate this approximation using the method of scalar 
variables, which in turn allows us to extract the general T dependence of most of the 
thermodynamic observables in an approximate form. Here, we apply this formalism 
to the computation of (/s. We explain the method of scalar variables in Appendix 
and the details of the adiabatic approximation in the scalar variables are given in 
Appendix y. 

For the scalar variable X (see Appendix ^ for a definition), the adiabatic ap- 
proximation means 

In Appendix |C| we present an independent argument based on the Einstein's equa- 
tions in scalar variables, for why this approximation holds in certain regimes. The 



-'^^Various coefficients in tliese equations differ from due to our different normalization of the 
dilaton kinetic term. 
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Figure 6: Comparison of the exact C,/ s with the adiabatic approximation in the variable 
\h- Solid(red) curve is the full numerical result and the dashed(blue) curve follows from 
( 13:231) . 



fluctuation equation ( p.6| ) greatly simplifles with ( |3.22| ). In fact, as shown in Ap- 
pendix y, the solution becomes independent of (f). With unit normalization on bound- 
ary, the adiabatic solution in the entire range of G (— cxd, (p^) becomes hadb{4>) = 1- 



Consequently, the coefficient Ch in (|3.12|) becomes unity, hence: 



3 ( V'M 
3271 V Vi^h) 



(3.23) 



We plot this function in Xh in figure ^, where we also provide the exact numerical 
result for comparison. Note that in figure ^ the whole large black-hole branch has 
been compressed at the left of the figure for < 0.04 The same functions in the 
variable T/Tc are plotted in figure |^. 

The validity of the adiabatic approximation equation ( 3.22 ), is determined by 
the rate which V /V varies with 0. In particular, the approximation becomes exact 
in the limits where V /V becomes constant. This happens for a constant potential 
or a potential that is a single power of A (eponential in (p). This is the case in the 
UV (0 —7- — oo, where the potential becomes a constant) and the IR (0 — t- +oo where 
the potential becomes a power law.) . Therefore equation (|3.23|) allows us to extract 
the analytic behavior of C/s in the limits (ph — )■ ±oo. 

The numerical values one obtains from ( p.23| ) in the intermediate region may 
differ from the exact result ( |3.12| ) considerably, especially near T^. However, we 
expect that the general shape will be similar. We refer to Appendix ^ for further 
details. 

Finally, the adiabatic approximation hints at why, in the particular background 
that we study, (/s at Tc is close to the limit value ( ^.16 ): In order to see this we 




Figure 7: Comparison of tlie exact C,/s with the adiabatic approximation in variable T. 
Solid(blue) curve is the full numerical result and the dashed(red) curve follows from ( 3.23| ). 



rewrite ( p.23| ) as 

c 

s 



—X\ (3.24) 

37r 



adb 

In the limit (|3.16|) we have X — )■ —1/2. The only other point where X = 1/2, is at 
the minimum of the string frame scale factor 0^,. This is the point where the confining 
string saturates ^9] . On the other hand, we expect on general physical grounds that 



the de- confinement phase transition happens near this point, i.e. 0^ ~ 0*- Thus, the 
adiabatic formula predicts that C/^i.'t'c) be close to the limit value l/Gvr.^'^ 

3.4 Buchel's bound 

In [^, Buchel proposed a bound for the ratio of the bulk and shear viscosities, moti- 
vated by certain well-understood holographic examples. In 4 space-time dimensions 
the Buchel bound reads. 



l>2^--cij, (3.25) 

We note that the bound is proposed to hold in the entire range of temperature from 
Tc to oo. This bound is trivially satisfied for exact conformal theories such as A/" = 4 



YM, and saturated in theories on Dp branes With the numerical evaluation 

at hand, we can check ( |3.25| ) in our case. In figure § (a) we plot the LHS and RHS 
of the bound We clearly see that the bound is satisfied for all temperatures. As 



"'^'^This argument may break down for two (dependent) reasons: First of all the adiabatic approx- 
imation becomes lees good near (pc- This is because, in this region V' /V varies relatively more 
rapidly as a function of (j). Secondly, precisely because of this, even though (j)c is not far away from 
(/)* the difference can result in a considerable change in the value of C/s through ( |3.23 ). 
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Since this theory contains two derivatives only, - has the universal value I/A-k. 
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Figure 8: (a) Comparison of C,/ri (solid line) and the RHS of ( |3.25| ) (dashed line), obtained 
using the speed of sound of the IHQCD model [^]. (b) Plot of the function C(T) defined in 
equation (3.26) as a function of temperature. The horizontal dashed line indicates where 
Buchel's bound is saturated. We see that the bound is satisfied in the entire range of 
temperatures. 



expected, both the LHS and the RHS of ( |3.25| ) vanishes in the high T conformal 
limit. 

A clear picture of Buchel's bound is obtained by defining the function: 



C{T) 



2(l/3-c2)^ 



(3.26) 



in terms of which the bound is simply C > 1. In Figure ^ (b) we show the function 
C(T) obtained numerically in our IHQCD model, between T^. and 5Tc. The values 
of this function are mildly dependent on temperature, and are between 1.5 and 2, 
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the same range of values that were recently considered in the hydrodynamic codes 
by Heinz and Song |^ . 



4. The drag force on strings and heavy quarks 

We will now consider an (external) heavy quark moving through an infinite volume 
of gluon plasma with a fixed velocity f at a finite temperature T Q . The quark 
feels a drag force coming from its interaction with the plasma and an external force 
has to be applied in order for it to keep a constant velocity. In a more realistic set 
up one would like to describe the deceleration caused by the drag. 

The heavy external quark can be described by a string whose endpoint is at the 
boundary. One can accommodate flavor by introducing D-branes, but we will not do 
this here. A first step is to describe the classical string "trailing" the quark. 

We consider the Nambu-Goto action on the worldsheet of the string. 

Sng = J dadr^det {-gMnd^X^ dpX^) , (4.1) 

where the metric is the string frame metric. The ansatz we are going to use to 



describe the trailing string is, p2|, 



X^=vt + i{r), X2 = X^ = 0, (4.2) 
along with the gauge choice 

a = r, r = t , (4.3) 

where r is the (radial) holographic coordinate. The string is moving in the X^ 
direction. 

This is a "steady-state" description of the moving quark as acceleration and 
deceleration are not taken into account. For a generic background the action of the 
string becomes 

S = j dtdr ^y-googrr - 9oo9iiC'^ - giigrrv"^ ■ (4.4) 

Note that g^Q is negative, and we should check whether our solution produces a real 
action. For example a straight string stretching from the quark to the horizon is a 
solution to the equations of motion but has imaginary action. 

We note that the action does not depend on ^ but only its derivative, therefore 
the corresponding "momentum" is conserved 



We solve for ^' to obtain 



V googu (1 + googu/ i2n£1ni:y) 



(4.6) 



The numerator changes sign at some finite value of the fifth coordinate r^. For the 
solution to be real, the denominator has to change sign at the same point. We 
therefore determine via the equation 



goo{rs) + gii{rs)v'^ = , 



and the constant momentum 

2 ^ goo{rs)giiirs 
^« - (27r£^)2 

Writing the string-frame metric as 



(4.7) 



(4. 



2 „2As 



ds = e 



T 



— / dt^ + dx ■ dx 



( ^y?!) becomes 



= f{rs) 

The induced world-sheet metric is therefore 



(4.9) 



(4.10) 



yap ^ 1 „2As(rs),,2 



f{r) Y e4'4s{'-)/(r)-e4'4i,{rs)„2 

g4As(r) j2(-j,-)_^4g4As{rs) 



/(r) Y e4'4s('-)/(r)_e4'4s('-s)«2 J2 j-^) |^g4As (r) j(-^)_^2g4As(r3) j 



(4.11) 



We can change the time coordinate to obtain a diagonal induced metric t 
T + C{r) with 

g2As (rs)^2 



fir)^/{fir) - t;2)(e4^»M/(r) - e^^^'^h'^) 



The new metric is 



ds^ = e^^^W 



-(/(r)-t;^)rfr^ + 



(e4^«W/(r) - e4^''('"«)t;2) 



(4.12) 



and near r = it has the expansion 



32A4rs) 



(4t;2A;(r,) + /'(r,))(r-r, 



+ 0(1) 



(4.13) 



This is a world-sheet black-hole metric with horizon at the turning point r = 



- 25 - 



4.1 The drag force 



The drag force on the quark can be determined by calculating the momentum that 
is lost by flowing along the string into the horizon: 



drag 



dpi ^ 1 9oo9u^' 

dt 27r^2 



(4.14) 



This can be obtained by considering the world-sheet Noether currents n^j'[Q and 
expressing the loss of momentum as AP^^ = / Hp This may be evaluated at any 
value of r, but it is more convenient to evaluate it at r = Ts. 
We flnally flnd that 



drag 



27r£2 



-9oo{rs)gu{rs) ■ 



(4.15) 



Using the form 
finally obtain 



of our finite-temperature metric in the string frame we 



o'2As{rs) 



drag 



V7(^ 



27r£2 



4/3 



(4.16) 



where in the second equality we expressed the force in terms of the Einstein-frame 
scale factor and the "running" dilaton. Substituting from ( [l.l(j| ) we obtain 



drag 



27r£2 



27r£2 



(4.17) 



Before proceeding further, we will evaluate the drag force for the conformal case of 
A/" = 4 SYM where 



= firs) = 1 - (vrTr,)^ 



(4.18) 



Substituting in ( [4.17] ) we obtain, pT|- 



conf 



(4.19) 



Moving on to YM, to compute the drag force from equation ( ]4.17| ) we must first 
determine ig in the IHQCD model. In this setup there is no analog of the Af = 4 
SYM relation ( 4.18| ) between i, ig and A. Rather, the fundamental string length ig 
is determined in a bottom-up fashion, by matching the effective string tension to the 
QCD string tension cTc derived from the lattice calculations. From (|2.15|) 



n- - _J_p2AsAr*) - _J_p-iMr,)\ ^4/3 



(4.20) 
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where is the point where the zero-temperature string scale factor As^o{t) has a 
minimum. For a typical value of (Tc ~ (440 MeVY 0] we find 



0.15 



(4.21) 



where £ is the radius of the asymptotic AdS space. 

On the other hand, unlike in A/" = 4 SYM, in the IHQCD model the value of the 
coupling A(rs) in equation (|4.17|) is not an extra parameter to be fixed by hand, but 
rather it is determined dynamically together with the background metric. 
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Figure 9: In this figure the ratio of the drag force in improved holographic QCD to the 
drag force in = 4 SYM is shown. The ratio is computed for different velocities as a 
function of temperature. The 't Hooft coupHng for the A" = 4 SYM theory is taken to be 
5.5. We chose this value as it is considered in the central region of possible values for the 
't Hooft coupling. It is seen that as the velocity increases the value of the ratio decreases. 
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Figure 10: In this figure the ratio of the drag force in improved holographic QCD to the 
drag force in A" = 4 SYM is shown. The ratio is computed for different temperatures as a 
function of velocity. The 't Hooft coupling for the A" = 4 SYM theory is taken to be 5.5. 
As temperature increases the value of the ratio decreases. 
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4.2 The relativistic asymptotics 

When f — 7- 1 then — )■ and we approach the boundary. Near the boundary (r — )■ 0) 
we have the foUowing asymptotics of the scale factor and the 't Hooft couphng, [H5| 



fir) 
and 



l + O 



A(r) 



1 



log(Ar^ 



1 



.Mr) 



O 



1 



log(Ar 
(4.22 



+• 



/3o log(rA) 
where is the position of the horizon. 

We therefore obtain for the turning point 



C(log(rA) 



-2\ 



(4.23) 



and the drag force 



l + O 



log(l — f 2) 



A(r, 



/So log 



+ ■ ■ 



drag 



We also use 



27r£2 v/r^ 
s(T) 457r£3s(T) 



+ 



47rM3 m 



(4.24) 



(4.25) 



(4.26) 



where s{T) the entropy per unit three- volume, and we write the relativistic asymp- 
totics of the drag force as, 



drag 



(-flog [i-v^]y 



- + 



(4.27) 



'45 Ts{T) 



Am 



- + ■ 



The force is proportional to the relativistic momentum combination v / \/l — v"^ 
modulo a power of log [1 — v"^]. This factor is present because, as argued in the 
asymptotic metric is AdS in the Einstein frame instead of the string frame. Its effects 
are not important phenomenologically. We discuss this issue further in Appendix ^ 

4.3 The non-relativistic asymptotics 

We now consider the opposite limit, f — 0. In this case the turning point asymptotes 
to the horizon, Vg — )■ Vh and we have the expansion 

o,2 



/(r) ~ 4TTT{rh - r) + 0{{rh - r) ) , Vs = Vh 



AtxT 



(4.28) 
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and 



(4.29) 



where primes are derivatives with respect to the conformal coordinate r. 
4.4 The diffusion time 

For a heavy quark with mass Mg we may rewrite ( [4.19[ ) as 



^ dp 1 MgV 

Fconf^-^ = —p , p=^=^ (4.30) 
dt T V 1 — V 

where the first equation defines the diffusion time r. In the conformal case, the 
diffusion time is constant, 

This is not anymore the case in QCD, where r defined as above is momentum de- 
pendent. We may still define it as in ( |4.3CI| ) in which case we obtain the following 
limits 



4 



£2 / A^2 \l 27r 
lim T = Mg^{ + ■■■ (4.33) 

4.5 Including the medium-induced correction to the quark mass 

In order to estimate the diffusion time of a quark of finite rest mass, we must take into 
account the fact that the mass of the quark receives medium- induced corrections. In 
other words, the mass appearing in equation ([4.30|) is a temperature-dependent quan- 
tity, Mg{T) 7^ Mg{T = 0). The ratio Mg{T)/Mg{0) can be estimated holographically 
by representing a static quark of finite mass by a static, straight string stretched 
along the radial direction starting at a point r = 7^ 0. At zero temperature, the IR 
endpoint of the string can be taken as the "confinement" radius, r^,, where the string 
frame metric reaches its minimum value; At finite temperature, the string ends in 
the IR at the BH horizon^^. The masses of the quark at zero and finite T are related 
to the worldsheet action evaluated on the static solution (r = t, cr = r) : 



15 



This representation ignores the fact that the kinetic mass of a moving quark may be different 



from the static mass [21 . We plan to treat this in the future. 

^^It would stop at the confinement radius if the latter were closer to the boundary than the 
horizon, i.e. if ?'*(T) < rh{T). However, in the model we are considering, in the big BH branch we 
find that the relation < is always satisfied. 
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Figure 11: In this figure the ratio of the diffusion time in the Improved Holographic QCD 
model to the diffusion time in = 4 SYM is shown. The 't Hooft coupling for = 4 
SYM is taken to be A = 5.5. The heavy quark has a mass of Mq = 1.3GeV. Note that 
with the definition of the diffusion time in ( [4.30| ) the ratio is the inverse of the ratio of the 
forces. A similar plot is valid for the bottom quark as well, as the mass drops out of the 
ratio, although the energy scales are different. In this plot the x-axis is taken to be in MeV 
units. As temperature increases the ratio also increases. 



= 2^ i rfre2^°WAy^(r) , M,(T) = ^ dr e'^^^^\'l\r) . (4.34) 

The value can be fixed numerically by matching Mg(0) to the physical quark 
mass, and translating the fundamental string tension in physical units by using the 
relation { W^ , with = {UOMeVf. This makes Mg{T) a function of Mg(0). The 
ratios Mg{T)/Mg{0) we found numerically in the model under consideration is shown 
in figure for the Charm (M(0) = 1.5GeV) and Bottom (M(0) = A.^GeV) quarks. 
The fact that, in the deconfined plasma, the quark mass decreases with increasing 
temperature is a direct consequence of the holographic framework^^, since for higher 
temperature, the distance to the horizon is smaller. An indication that this result 
may be in the right direction comes from the lattice computation of the shift in 
the position of the quarkonium resonance peak at finite temperature [Q: in the 
deconfined phase the charmonium peak moves to lower mass at higher temperature. 
Our result for the medium-induced shift in the constituent quark mass is consistent 
with these observations. 

We can now write the diffusion time from eqs. (|4.17|) and ( [4.30[) as: 



(Tca/I - v"^ \ A(r 



17 



For a possible field theoretical explanation of this phenomenon, see 
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M(T) 




Figure 12: Ratios between the thermal mass and the rest mass of the Charm (cm've 
labeled "c") and Bottom (curve labeled "b" ) quarks, as a function of temperature. 



where once again we have ehminated the fundamental string length using equation 
( [4.20|) . Given a set of zero- and finite-temperature solutions, equation (|4.35|) can be 
evaluated numerically for different values of the velocity and different quark masses. 
The results for the Charm (Mg(0) = 1.5 GeV) and Bottom (M = 4.5 GeV) quarks 
are displayed in figure |T3|. 
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Figure 13: Diffusion time for the Charm and Bottom quarks, as a function of energy, for 
different ratios of the temperature to the IHQCD transition temperature Tc- 



4.6 Temperature matching and diffusion time estimates 

An important question is how we should choose the temperature in our holographic 
model in order to compare our results with heavy-ion collision experiments. This is 
nontrivial, since our setup is designed to describe pure SU{Nc) Yang-Mills, whereas 
at RHIC temperatures there are 3 light quark flavors that become relevant. As a 
consequence, the critical temperatures and the number of degrees of freedom of the 
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two theories are not the same: for pure SU{Nc) Yang Mills we have A^^ — 1 degrees 
of freedom and a critical temperature around 260 MeV; For SU{Nc) QCD with 
Nf flavors the number of degrees of freedom is — 1 + N^Nf, and the transition 
temperature is lower, around 180 MeV. In our holographic model, the transition 
temperature in physical units was estimated to be Tc = 247 MeV P9[, i.e. close to 



the lattice result for the pure YM deconfining temperature. From now on, this is 
the value we will mean when we refer to T^. This is also close to the temperature of 
QGP at RHIC, which we will denote Tqgp, and is estimated to be around 250 MeV. 
Since this value is uncertain, below we give our results for a range of temperatures 
between 200 MeV and 400 MeV. The higher temperatures will be relevant for the 
LHC ion collision experiments (see e.g. 16^). 

Based on these considerations, there are different ways of fixing the temperature 
(see e.g. the recent review [l^): one direct and two alternative schemes (that we call 
the energy and entropy scheme). 

• Direct scheme: The temperature of the holographic model is identified with 
the temperature of the QGP in the experimental situation (at RHIC or LHC), 

rp{dir) rp 

J-ihqcd - J-QGP- 

• Energy scheme: One matches the energy densities, rather than the tempera- 
tures. The energy density at RHIC is approximately (treating the QCD plasma 
as a free gas^^.) egap ^ (vrVl5)(iV2 - 1 + N,Nf){TQGpY. For N, = Nf = 3, 
asking that our energy density matches this value requires us to consider the 
holographic model at temperature T^f} , given by 



^^hUTlLa) ^ 11-2(Tqgp)' (4.36) 



• Entropy scheme: Instead of matching the energy densities, alternatively one 
can match the entropy density s, which for the QGP, in the free gas approxi- 
mation, is given by aqcp ~ A'n'^/A'b{Nl - 1 + NcNf){TQGpY. This leads to the 
identification: 

s^Hca{TiulJ = 14.9(Tqgp)' (4.37) 

The temperature translation table between the various schemes is shown in Table 
|I|. In that table, = 247Mel^ is the deconfining temperature of the holographic 
model. 

In Figure ^we show the comparison between the diffusion times, as a function of 
initial quark momentum, in the different schemes for the Charm and Bottom quarks, 
at the temperature Tqgp = 250MeV. 

"'^^This is itself an approximation, since as we know both from experiment and in our holographic 
model, the plasma is strongly coupled up to temperatures of a few T^, 
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Tqgp (MeV) 


Tqgp/Tc 


T]^lcA (MeV) 


rp{t) Irp 

ihqcdl -'- c 


TlCcd (MeV) 


ihqcdl c 


190 


0.77 


259 


1.05 


274 


1.11 


220 


0.89 


290 


1.18 


302 


1.23 


250 


1.01 


325 


1.31 


335 


1.35 


280 


1.13 


361 


1.46 


368 


1.49 


310 


1.26 


398 


1.61 


402 


1.63 


340 


1.38 


434 


1.76 


437 


1.77 


370 


1.50 


471 


1.90 


472 


1.91 


400 


1.62 


508 


2.06 


507 


2.05 



Table 1: Translation table between different temperature identification schemes. The first 
two columns display temperatures in the direct scheme, (in which the temperature of the 
holographic model matches the physical QGP temperature) and the corresponding ratio 
to the IHQCD critical temperature, that was fixed by YM lattice results at Tc = 247 MeV 
1 49]. The third and fourth columns display the corresponding temperatures (and respective 
ratios to Tc) in the energy scheme, and the last two in the entropy scheme. 




Figure 14: Diffusion times for the Charm and Bottom quarks, as a function of initial 
momentum, at Tqgp = 250 MeV. The different lines represent the in the direct scheme 
(solid), energy scheme (dashed) and entropy scheme (dash-dotted), all corresponding to 
the same temperature Tqgp = 250 MeV. 



The results for the diffusion times at different temperatures, computed at a 
reference heavy quark initial momentum p \Q GeV, are displayed in Tables ^ and 
^. We see that there is little practical difference between the entropy and energy 
schemes; on the other hand the difference between the direct scheme and the two 
alternative schemes can be quite substantial. 
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TQGP,MeV 


Tcharm (fm/c) 

(direct) 


Tcharm (fm/c) 

(energy) 


r^harm (fm/ C ) 

(entropy) 


220 


- 


4.0 


3.6 


250 


5.7 


3.1 


3.0 


280 


4.3 


2.6 


2.5 


310 


3.5 


2.1 


2.1 


340 


2.9 


1.8 


1.8 


370 


2.5 


1.5 


1.5 


400 


2.1 


1.3 


1.3 



Table 2: The diffusion times for the charm quark are shown for different temperatures, 
in the three different schemes. Diffusion times have been evaluated with a quark initial 
momentum fixed at p w 10 GeV. 



TQGp{MeV) 


Tbottom (fm/c) 


Tbottom (fm/c) 


Tbottom (fm/c) 




(direct) 


(energy) 


(entropy) 


220 




8.9 


8.4 


250 


11.4 


7.5 


7.1 


280 


10.1 


6.3 


6.1 


310 


8.6 


5.4 


5.3 


340 


7.5 


4.7 


4.7 


370 


6.6 


4.1 


4.1 


400 


5.8 


3.6 


3.6 



Table 3: Diffusion times for the bottom quark are shown for different temperatures, 
in the three different schemes. Diffusion times have been evaluated with a quark initial 
momentum fixed at p ~ 10 GeV. 



5. Jet quenching parameter 

In this Section we discuss the jet quenching parameter in the class of holographic 
models under consideration, and we estimate its numerical value for the concrete 
model with potential (|2.16| ) and parameters fixed as in For the holographic 



computation, we will follow |]T9|, 0. There is another method available |^8[, but we 
will not use it here. 

The jet-quenching parameter q provides a measure of the dissipation of the 
plasma and it has been associated to the behavior of a Wilson loop joining two 
light-like lines. We consider two light-like lines which extend for a distance L~ and 
are situated distance L apart in a transverse coordinate. Then q is given by the large 

behavior of the Wilson loop 

W ~ e"^''"^"^' . (5.1) 
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We consider the bulk string frame metric 



ds^ = e^^^^'-) (^-f{r)df + dx^ + . 



(5.2) 



To address the problem of the Wilson loop we make a change of coordinates to light 
cone coordinates for the boundary theory 



x'^ — Xi + t X — Xi — t 



(5.3) 



for which the metric becomes 



ds'^ = e^^» (dxj + dxl + ^{l- f){dxl + dxl) + (1 + f)dx+dx- + . (5.4) 

The Wilson loop in question stretches across X2, and lies at a constant x+,X3. It is 
convenient to choose a world-sheet gauge in which 

X- = X2 = o . (5.5) 

Then the action of the string stretching between the two hues is given by 

S-^J dadr^ -det{gMNdaX^dpX^) (5.6) 
and assuming a profile of r = r{a) we obtain 



I 



dx2 e 



2As 



(1-/) 



1 + 



„/2 



/ 



(5.7) 



The integrand does not depend explicitly on X2, so there is a conserved quantity, c: 

,dC „ c 



^ dr' ^ 



V2 



which leads to 



(5.8) 



(5.9) 



A first assessment of this relation involves determining the zeros and the region of 
positivity of the right-hand side. / is always positive and vanishes at the horizon. 
For the second factor we need the asymptotics of e'^^'(l — /). This factor remains 
positive and bounded from below in the interior and up to the horizon. It vanishes 
however logarithmically near the boundary as 



e4^«(l-/) = 7rm3^('''') 



3o log(Ar) 



log(Ar) 



(5.10) 
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Figure 15: In this figure the combination (1 — f)e " is plotted as a function of the 
radial distance, for several temperatures. The radial distance is given in units of the 
horizon position for the black hole at the critical temperature T^. All curves stop at the 
corresponding horizon position. 



This is unlike the conformal case where we obtain a constant 

AAs 



e*^=(l-/) ={ixT(.f (5.11) 

conformal 



The behavior in (|5.1CI| ) is a model artifact and is analyzed in appendix ^ 



Because of this, for fixed c, there is a region near the boundary where r'^ becomes 
negative. At this stage we will avoid this region, by using a modified boundary at 
r = e. We will later show that this gymnastics will be irrelevant for the computation 
of the jet quenching parameter, as it involves effectively the limit c — )■ 0. 

We will place the modified boundary r = e a bit inward from the place r = r^m 
where the factor - — ^ — 1 vanishes: 



e'^'^'--\l-f{rrmn)) = c^ (5.12) 

Therefore we choose Tmin < e- 

Then, in the range e < r < Vh the factor - — -^3—^ — 1 is positive for sufficiently 
small c. In this same range, r' vanishes only at r = r^. This is the true turning point 
of the string world-sheet. This is also what happens in the conformal case. It is also 
intuitively obvious that the relevant Wilson loop must sample also the region near 
the horizon. 

The constant c is determined by the fact that the two light-like Wilson loops are 
a X2 = L distance apart. 

L _ cdr 

2 "A ^/(e4A.(i_/)_=^- ^^-^^^ 
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The denominator vanishes at the turning point. The singularity is integrable^^. 
Therefore, as we are interested in the small L region, it is obvious from the expression 
above that that c must also be small in the same limit. 
This relation can then be expanded in powers of c as 

+ W ' +0{c'). (5.14) 



2c Je /7(w) yjiwF 

Therefore to leading order in L 



2/. 



0{L^) (5.15) 



We are now ready to evaluate the Nambu-Goto action of the extremal configu- 
ration we have found. Starting from (|577|) , we substitute r' from ( ^.9|) , and change 
integration variable from X2 t to obtain 

^ '^-r - , (5,16) 

27,11], V2/(e"'(l -/)-<:') 

As in [19, 20 1, we subtract from equation (|5.16|) the action of two free string 
straight worldsheets that hang down to the horizon. To compute this action a con- 
venient choice of gauge is x_ = r, r = a. The action of each sheet is 



S„ = ^f *e-yi^ (5.17) 

The subtracted action is therefore: 

Sr = S-2So = —- + Ojc") , (5.18) 



Using now (|5.15|) to substitute c we finally obtain 



Sr = + 0{L^) . (5.19) 



So far we have evaluated the relevant Wilson loop in the fundamental represen- 
tation (by using probe quarks). On the other hand, the Wilson loop that defines the 
jet-quenching parameter is an adjoint one. We can obtain it in the large- A'^c limit from 
the fundamental using tr Adjoint = ^^Fundamental- ^6 finally extract the jet-quenching 
parameter as 

V2 1 



^ = h-.—. Tr • (5-20) 



-'^^Even if we choose e = r'mm) the new singularity at r = is also integrable as suggested from 
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We are now ready to remove the cutoff. As the integral appearing is now well-defined 
up to the real boundary r = we may rewrite it as 



e ^^^c/r 



e ^^'•dr 



lie) , lie) 



e ^^'•dr 




(5.2i; 



In appendix we obtain the small e estimate of /(e) that vanishes as ~ e(log e) 3 . 
We may finally write^° 

^ = ^— Kr • (5.22) 
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Figure 16: In this figure the ratio of the jet quenching parameter in our model to the 



jet quenching parameter in A/" = 4 is shown. The integral present in equation ( 5.20 ) has 
been numericahy calculated from an effective cutoff at r = r/j/1000. The jet quenching 
parameter in = 4 SYM has been calculated with X'tHooft = 5.5. 



From equation ( 5.22 ) we obtain, in the conformal case: 



Q'conformal 



(5.23) 



The conformal value, for the median value of A = 5.5 and T ^ 250 MeV gives 
•^conformal — 1.95 GcV^/fm where we used the conversion 1 GeV~ 5 fm~^. 



Numerical evaluation of equation ( ^.22|) in the non-conformal IHQCD setup^^ 
gives us a value of q which is lower (at a given temperature) than the conformal value. 



^'^In practise, the previous discussion including regularizing the UV is academic. The numerical 
calculation is done with a finite cutoff where the boundary conditions for the couplings are imposed. 



^^In this case, the value of £s appearing in equation (5.22) is fixed as explained in Section 4. 
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Figure 17: The jet quenching parameter q for the Improved Holographic QCD model and 
A/" = 4 SYM is shown in units of GeV'^ / fm for a region close to T = Tc. The smallest 
dashed curve is the ihQCD result with an effective cutoff of r cutoff = r/j/lOOO. The small 
dashed curve is the ihQCD result with the cutoff from the mass of the Bottom quark. The 
medium dashed curve has a cutoff coming from the Charm mass and and largest dashed 
curve is the = 4 SYM result. 





T/Tc 



Figure 18: The jet quenching parameter q for the Improved Holographic QCD model 
(lower curve) and A/" = 4 SYM (upper curve) are shown in units of GeV"^ / fm for temper- 
atures up to r = ATc- 



as shown in Figures To, T7 and la. Tables ^ to |^ display the numerical values of 
the jet quenching parameter at different temperatures in the experimentally relevant 
range, in different temperature matching schemes. 



6. Discussion and summary 

In this paper we have examined several aspects associated with the physics of ther- 
mal transport phenomena in gluon plasma with potential applications to heavy ion 
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TQGP,MeV 


q {GeV^/fm) qi {GeV^/fm) 
(direct) (direct) 


220 


- 


250 


0.5 0.6 


280 


0.8 0.8 


310 


1.1 1.1 


340 


1.4 1.4 


370 


1.8 1.8 


400 


2.2 2.2 



Table 4: This table shows the jet quenching parameter q computed with different cutoffs 
for the different temperatures shown in the first column. The computation is done in the 
direct scheme. The second column shows q with a cutoff at r cutoff = '''/i/lOOO, where r/j 
is the location of the horizon. In accordance with the conclusions of appendix q does not 
change significantly as we vary the cutoff from r/i/1000 to r/^/lOO. 



Tqgp, MeV 


q {GeVyfm) 


q {GeVyfm) 


q (GeV^/fm) 




(direct) 


(energy) 


(entropy) 


220 




0.9 


1.0 


250 


0.5 


1.2 


1.3 


280 


0.8 


1.6 


1.7 


310 


1.1 


2.1 


2.2 


340 


1.4 


2.7 


2.8 


370 


1.8 


3.4 


3.4 


400 


2.2 


4.2 


4.2 



Table 5: This table displays the jet quenching parameter q using the three different 
comparison schemes. For lower temperatures the "entropy scheme" gives higher values. 
As energy is increased the energy and entropy schemes temperatures start to coincide and 
there is little difference in the jet quenching parameter as well. 



collisions. We have used as basic model for our calculations the 5D Einstein-dilaton 
model with a potential proposed in This is a hybrid model that incorporates 



features coming from string theory as well as features originating in YM, [^]. We 
have also used a potential, whose two phenomenological parameters have been fit to 
lattice YM data, This model, named Improved Holographic QCD is expected 
to be a very good approximation to several aspects of YM physics. 

In this context we calculated the bulk viscosity by calculating the low-frequency 
asymptotics of the appropriate two-point function of the energy-momentum tensor. 
We have further calculated the drag force on a heavy quark by extending the dragging 
string calculation done earlier in the context of A/" = 4 SYM, ^ |2^. Finally 



we have calculated the jet-quenching parameter q defined by the expectation value 
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Tqgp, MeV 


qcharm (GeV^/fm) 


qcharm (GeV^/fm) 


qcharm (GeV'^/fm) 




(direct) 


(energy) 


(entropy) 


220 


- 


1.3 


1.5 


250 


0.8 


1.8 


2.0 


280 


1.2 


2.6 


2.8 


310 


1.7 


3.5 


3.6 


340 


2.2 


4.6 


4.7 


370 


2.8 


5.9 


6.0 


400 


3.6 


7.6 


7.5 



Table 6: This table displays the jet quenching parameter q using the three different 
comparison schemes with an effective cutoff provided by the mass of the Charm quark. 
Again, for lower temperatures the "entropy scheme" gives higher values. As energy is 
increased the energy and entropy schemes temperatures start to coincide and there is little 
difference in the jet quenching parameter as well. Also when the temperature approaches 
the quark mass the picture of the heavy quark as a hanging string collapses and results are 
not reliable. 



Tqgp, MeV 


qbottora {GeV'^/fm) 


qbottora (GeV'^/fm) 


qbottora (GeV'^/fm) 




(direct) 


(energy) 


(entropy) 


220 




1.0 


1.1 


250 


0.6 


1.4 


1.5 


280 


0.9 


1.9 


2.0 


310 


1.2 


2.5 


2.6 


340 


1.6 


3.2 


3.2 


370 


2.0 


4.0 


4.0 


400 


2.5 


5.0 


4.9 



Table 7: This table displays the jet quenching parameter q using the three different 
comparison schemes with an effective cutoff provided by the mass of the Bottom quark. 
The results are close to the q results computed in Table ^ since the mass of the Bottom 
quark is much larger than the temperatures we examine. 



of a light-like Wilson loop, by adapting the calculation of |19, BQ] from A/ = 4 



Unlike the case of A/" = 4, our calculations here are numerical as no analytical 
solutions are known for Improved Holographic QCD with the appropriate potential. 
We have however derived analytically various asymptotics of the results relevant for 
high energy, low velocity or high temperature. 

Before discussing the results, it is appropriate at this point to take a critical look 
and analyze potential sources of (systematic) error in our calculations. 

• Holographic models are reliable in the context of large-A^^c expansion of the 
SU{Nc) gauge theory. Therefore a priori, our results should be understood as 
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the leading order 0{1) part in the large-Nc expansion. 

The issue, however, is a bit more comphcated by the fact the the model we are 
employing is semi-phenomenological and therefore contains two phenomeno- 
logical parameters (apart from the ones expected in YM) that have been fit to 
data in [^. Although several lattice data are known at large-Nc, |6^, others 



are not. In particular, the detailed thermodynamics of large-Nc YM is currently 
being calculated on the lattice, Therefore not all relevant input data we 



used have been computed at large Nc- 

In this sense the semi-phenomenological model we are using is positioned some- 
where between Nc=3, YM and Nc = oo YM. It is known so far that the dif- 
ference in many observables in the gluon sector between these two points is of 
the order of 5% or less. 

There are no dynamical flavor degrees of freedom incorporated in the model 



used. In |38, 39 1 the incorporation of flavor branes was described at the semi- 
quantitative level. We have assumed that we work in the "quenched" approx- 
imation: the number of flavors Nf <^ Nc which implies in particular that 
fermions loops are suppressed by a factor of ^ 1. 

The configuration of flavor branes is expected to involve a pair of space-filling 
and per quark flavor. These branes enter at the AdS boundary and at 
some point in the interior they are expected to fuse signalling chiral symmetry 
breaking. The configuration in the broken phase involves a space-filling brane 
that folds on itself and resembles closely the branes described in |^ using 
Boundary CFT. 

The bare mass of the associated quarks enters as a source boundary condition 
on the relevant tachyon field, [^. The higher the mass the stronger is the 
tendency on the tachyon to diverge in the IR. In the deconfined phase we expect 
that branes associated to light quarks cross the BH horizon and this signals 
the melting of the associated mesons. Branes associated to heavy quarks, will 
fuse outside the horizon, signaling the stability of the associated mesons. These 
expectations are qualitative. They have been observed in toy models, ||6^ but 
have not been yet calculated in a reliable extension of the present setup. 



Our estimate is that the bare quark mass is related to the flavor brane position 
in the following way: The energy of a string stretching from to r^,, (the 
equilibrium position of string world-sheets) is equal to the bare quark mass, as 



detailed in section [4.5| . This is expected to be asymptotically correct when the 
mass of the quark is much larger than the dynamical scale of the gauge theory, 
and we therefore do not expect a large source of error for charm and bottom 
quarks. 
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There are several other sources of error, that enter between using the quantities 
computed here and comparing them to the eventual experimental data. Most 
of them have been described preciously, ^ Q , and we do not have 



much more to add here. We would like however to mention one extra important 
aspect: deciding the appropriate temperature to be used in comparisons with 
data. This is an issue because in YM the deconfining transition is first order 
(instead of the expected cross-over in the theory with quarks) with a transition 
temperature that is about 50% larger than in QCD. 

There is therefore a non-trivial comparison to be made. We do not know 
the best way to compare, but we have explored three different matchings: 
taking the same temperature, the same energy density or the same entropy 
density. Until a computation is made taking into account the fermionic degrees 
of freedom, this choice introduces an extra systematic error in the comparison. 

This ambiguity is the same one that arises when one fixes the temperature in 
the holographic computations using the A/" = 4 theory. In that case however, 
one must also fix the A/" = 4 coupling constant, and this introduces an extra 
source of error. In our model this is not an issue, since the coupling constant 
runs. All observable quantities we compute are independent of the value of the 
coupling at a given energy, thus we do not need to fix an extra parameter by 
hand. 

There are further limitations on the range of applicability of the drag force 



and jet-quenching calculations, that have been discussed in the literature ^0 



71[1 . In the drag-force calculation the velocity is limited by the position of the 
associated flavor brane. In a sense the world-sheet horizon should be kept away 
from that brane so that standard calculations of the drag force are reliable. The 
jet-quenching parameter seems valid in the opposite regime. 

There is a further important issue concerning the physics of heavy quarks in the 
context of the QGP. It has been argued from various points of view,[^],p4[-^l| that 



the motion of a heavy quark propagating and interacting in QGP is very similar to 
Brownian motion. The associated description starts with a Langevin equation which 
contains two ingredients: a classical force (the drag force) and a fluctuation force 
(the "noise") characterized by a a diffusion coefficient in the (late-time) Gaussian 
case. 

The distribution for the kinematic data then describes a Fokker-Planck equa- 
tion. In the standard non-relativistic case, the assumption of a Maxwell equilibrium 
solution to the Fokker-Planck equation provides a relation between the classical force 
and the stochastic force known in the simplest cases as the Einstein relation. 

One of the relevant ingredients in the case of QGP is that the description of the 
Brownian motion must be relativistic. Relativistic Langevin evolutions have been 



-43- 



described already in the relevant mathematical literature, |32|. However, in the early 
literature, it was assumed that the relativistic Maxwell distribution is an equilibrium 
solution to the Fokker-Planck equations. This leads to an Einstein-like relation 
that is problematic at high temperatures. This is taken as a hint that the initial 
assumption is false. Moreover, in the case of heavy-quark diffusion the longitudinal 
and transverse directions behave very differently. A recent series of papers, |p^8||-PT| 
derived the Langevin-type evolution of a heavy quark in the context of AdS/CFT 
by studying the small oscillations of the trailing string solution that describes the 
average motion of the quark. In particular, it was shown in that the fluctuating 



force is strongly influenced by the existence of an induced 2d black-hole metric and an 
associated world-sheet horizon in the semiclassical trailing string solution, as detailed 
in section^. In the non-relativistic limit, this world-sheet horizon and the bulk black- 
hole horizon coincide. The relevant diffusion coefficients are therefore computed from 
thermal two-point functions for the string fluctuations. 

The analogous computation of such thermal correlators in our case, is more in- 
volved than the JV = 4 case and will be reported in a future publication. 

Below we give a summary of our results. 

Bulk Viscosity: We have computed the bulk viscosity by calculating the low fre- 
quency asymptotics of the appropriate stress tensor correlator holographically. We 
find that the bulk viscosity rises near the phase transition but stays always below 
the shear viscosity. It floats somewhat above the Buchel bound, with a coefficient 
of proportionality varying between 1 and 2. Therefore it is expected to affect the 



elliptic flow at the small percentage level |p9| , |T3[. Knowledge of the bulk viscosity 
is important in extracting the shear viscosity from the data. This result is not in 
agreement with the lattice result near Tc. In particular the lattice result gives a value 
for the viscosity that is ten times larger. 

The bulk viscosity keeps increasing in the black-hole branch below the transition 
point until the large BH turns into the small BH at a temperature T^m- The bulk 
viscosity on the small BH background is always larger than the respective one in the 
large BH background. In particular, we showed that the T derivative of the quantity 
(/s diverges at T^m- This is the holographic reason for the presence of a peak in 



(/s near Tc. On the other hand, as it is shown in presence of T^m (i-e. a small 
BH branch) is in one-to-one correspondence with color confinement at zero T. We 
arrive thus at the suggestion that in a (large N) gauge theory that confines color at 
zero T, there shall be a rise in (/s near Tc. 

An important ingredient here was the value of the viscosity asymptotically in 
the small BH branch. There we correlated precisely its asymptotic value to the IR 
behavior of the potential. Taking also account the fact that this asymptotic value 
is very close to the value of the bulk viscosity near Tc, we can derive bounds that 
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suggest that the bulk viscosity cannot increase a lot near Tc. 

Drag Force: The drag force we have calculated has the expected behavior. Al- 
though it increases with temperature, it does so slower than in A/" = 4 SYM, signaling 
the effects of asymptotic freedom. This feature is at odds with the robustness obser- 
vations for the drag force of []72 . 



Diffusion Time: Based on the drag force calculation we have computed the dif- 
fusion times for a heavy external quark. The numerical values we obtain are in 



agreement with phenomenological models ^Tj. To accommodate for the fact that 
our models exhibits a phase transition around T = 247 MeV (i.e. about 30% higher 
than in QCD), we compare our results using alternative schemes, as proposed in [Q. 
For example, for an external Charm quark of momentum p = 10 GeV we find (in the 
alternative scheme) a diffusion time of r = 2.6 fm at temperature T = 280 MeV. 
Similarly, for a Bottom quark of the same momentum and at the same temperature 
we find r = 6.3 fm. Generally the numbers we obtain are close to those obtained 
by [O and m. 



Jet Quenching: We have also calculated the jet quenching parameter of this 
model, based on the formalism of [|19[ Q by computing the appropriate light-like 
Wilson loop. We find that q grows with temperature, but slower than the growth 
of J\f = 4 SYM result. Again this can be attributed to the incorporation of asymp- 
totic freedom in our model. Using the alternative scheme to compare with experiment 
we find that our results are close to the lower quoted values of q. For example, for a 
temperature of T = 290 MeV, which in the alternative "energy scheme" corresponds 
to a temperature of T = 395 MeV in our model, we find that q ^ 2GeV'^/ fm. 

However, the numbers obtained for this particular definition of jet quenching 
parameter seem rather low and indicate that this may not be the most appropriate 
definition in the holographic context. There are other ways to define q, in particular 
using the fluctuations of the trailing string solution. This is gives a direct and more 
detailed input in the associated Langevin dynamics and captures the asymmetry 
between longitudinal and transverse fluctuations. It would be interesting to compute 



this, along the lines set in [2S, 30, 31 and we are currently pursuing that aim. 



Note added in proof 

Since this paper has appeared in the archive, two papers appeared that have a direct 
connection to some of the issues discussed here. In reference [|7S| a high precision 
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lattice calculation of the thermodynamics was performed at various Nc- The results 
suggest that the thermodynamic functions vary very little with Nc although the phase 
transition becomes sharper as increases. The thermodynamic functions and in 
particular the trace-anomaly calculated from the Improved Holographic QCD model 
match very well the lattice data. 

In reference a detailed study of the hydrodymamics with a high-viscosity 
regime was performed. It was found that cavitation ensues for bulk or shear viscosity 



values a few times the PSS value, thus corroborating earlier numerical evidence, 



Our results indicate that cavitation (and therefore breakdown of the hydrodynamic 
description) is not expected to happen in the deconfined phase of the quark gluon 
plasma. 
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APPENDIX 



A. The scalar variables and evaluation of the bulk viscosity 



To determine ( we need to solve the fluctuation equation ( p^ numerically. This 
requires knowledge of the background functions A, B and / as functions of 0. A 
very convenient reformulation of the Einstein's equations, especially when the radial 
variable is taken as is explained in section 7 of that we review here. 



One can reduce the number of Einstein's equations by introducing the following 
scalar variables: 

XU) = = (A.l) 

where we defined / = exp((7). Note that X and Y are invariant under radial coordi- 
nate transformations. These variables obey the following first order equations: 

dX 4^^ ^2 , , 3 dlogV 



dY 4 Y 



As shown in , the thermodynamics of the dual field theory are completely deter- 
mined by knowledge of X and F as a function of 0. Roughly speaking, Y is dual to 
the enthalpy and X to the energy of the gluon fiuid. 



In solving ( [A.2| ) and (|A.3| ) one imposes the boundary conditions at the horizon. 



The regularity of horizon requires 



3 

as — )■ Solving ( |A.2|) near the horizon determines 



(ph- (p 

X ^ -\Yh + 0{(t>h-(t^), (A.4) 



-HI- 

Having solved for Y and X, one determines the metric functions A and / as, 

^ = ^o+f^d~^, (A.6) 
AY ~ 

g = \ogf = j (A.7) 

Now, let us compute the last metric function B. The metric written in the r- frame 
and the 0-frame are: 

ds^ = e'^'i-fdt' + dx' + Y^ = e'^i-fdt' + dx') + e^^^. (A.8) 
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Comparison determines, 



B = A-\og\'^\. (A.9) 



dr 

In the formulation of the scalar variables, d(j)/dr is given by. 



Thus, one finds B as 



= -e^ 3 J-^^. (A.io) 

dr I 



B = ^ I X -log|3X|. (A.ll) 



Having found the metric functions in X and Y variables, one can rewrite the 
fiuctuation equation (|3.6|) . There are various cancellations most notably in rewriting 
the w-dependent term in (|3.6| ): The temperature T is determined by the following 
equation, (see equation (H.67) pf [|^) in the scalar variables: 

y^n0o)e^o-4''#i_ (A. 12) 



Now, using the equations ( [A. 3D , ( |A.6D , ( [A.7| ) and ( |A.12| ), the oo dependent term can 
be simplified as 

\37iTX J 

With similar simplifications, the entire ( ^.6[ ) equation can be written only in terms 
of X and Y functions: 

h'l^ = c{(t))h[^ + d{(t))hn, (A. 13) 

where 

c(0) = ( - + — — ) , (A.14) 

2 



X V 3 2X V 



d{<P) = -1^(1 -X^ + Y){1 + -^^) - ( e-/:H. (A.15) 



To summarize: Given one computes the functions X and Y from (|A.2|) and (|A.3| ) 



and the temperature from ( [A.12|) . Given these data, one solves (|A.13|) numerically 



(with the boundary conditions explained below ( |3.6| )). 

In passing, we note that the equation ( |A.13 ) can be put in a Riccati form by 



the change of variables hn = exp(J h): h' + h? = ch + d whose general solution can 
be found iff one knows a special solution. It is presumably possible to find a special 
solution for simple potentials V . 
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B. Bulk viscosity in the limit of vanishing black-hole 



Here, we fill in the details of the computation that leads to equation (|3.16| ). This 
follows from ( |3.12| ) in the high T limit on the small BH {Xh — ?■ oo). We first show 
that, the fluctuation coefficient \cb\ goes to 1 in this limit. Cb is given by the value 
of hii{Xh) that follows from solving ( |A.13D with u = 0, and the boundary condition 

/lll(-CX)) = 1. 

In , it was shown that in the Xh ^ oo limit, the functions X and Y simplify. In 



particular X(A) — )■ Xo(A) where Xq corresponds to the zero T solution and Y{X) — )■ 
everywhere except X = Xh- In fact, one can show that Y is proportional to a delta 
function 6{X — Xh) in the limit Xh — )■ oo. Thus, from (|A.15|) we observe that d{(j)) 
vanishes in this limit for all values of X < Xh- In fact it also vanishes at Xh because 
the term 1 + 3V' /8XV vanishes as A = A/j — ?■ oo. Therefore the fluctuation equation 
simplifies to 

The solution with the aforementioned boundary condition is, 

hu{^) = 1 + C r dt e^--"°(*\ (B.2) 

J —oo 

The integration constant C is determined by the second boundary condition h'{(j)h) = 
0. On the other hand, c(0) is positive definite in the limit — )■ oo. This is because 
V /V approaches to 4/3 whereas Xq approaches to —1/2. Hence the only way to 
obey the condition is to set C = 0, hence hn = 1 for all values of A in the limit 
A/j — 7- oo. We checked that this is indeed the case by numerical analysis. 



C. The adiabatic approximation in scalar variables 



The approximate solution explained in ^6|] is given by eqs. ( p.20| ) and (|3.21|) . Us- 
ing logs oc A and (|A.6|) , we observe starting from ( p.2(]| ) that the approximation 
translated in scalar variables implies. 



3 ^(0) 
8 V{^) ■ 



(C.l) 



To verify that the second equation ( p.21 ) leads to the same conclusion, we may use 
|: 

log s - 3 log T oc -4 Z"^' X - 3 log V{(ph) ■ (C.2) 



equation (7.38) of [|5 



On the other hand, \ogV{(ph) oc -| /'^'^ X [gSf. Therefore, we verify that (|3]2T|) also 
leads to (Q. 
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Figure 19: Comparison of tlie scalar function X, the adiabatic approximation Xadb ^-nd 
the corresponding zero-T variable Xq. A/i = 1 is chosen to be 1. Solid (black), short- 
dashed(red) and the long-dashed (blue) curves correspond to the full numerical result X, 
the adiabatic approximation Xadb and the zero-T result Xq respectively. 



We compare both sides of equation (|C.1|) in figure |T9| for a large enough (so 
that a wider range can be compared). On this figure we also plot Xq (the variable 
X for the zero-T theory) for comparison. 

We will now proceed to understand the approximate formula ( |C.1|) indepen- 
dently. Suppose that V /V is a slowly varying function of (p. Then, we claim that 
we can write 

X(0) = X,,,(0) + 5(0), (C.3) 

where 5(0) is small w.r.t Xadb everywhere (this also means that 5'(0) is small every- 
where) . 

Substitution of in (Q gives^^, 



Kdb oc 5. (C.4) 

Therefore, the condition X ^ Xadb is equivalent to V /V is slowly varying with 0, 
namely the condition for the adiabatic approximation. This argument also shows 
that in the limits where V /V becomes constant, in particular near the boundary 
— 7- — oo, and near the singularity — +oo (for 0/j — t- oo), the approximation 
becomes exact. Figure ^ supports our arguments above by numerical evidence. 
Here we plot the ratio {(/s{exact) — (/s{adb)) /{(/ s) (exact), namely the difference 
between the exact (numerical) result and the adiabatic approximation (normalized 
by the exact value) and the function |(V/1^)'(0)|. The latter provides the criterion 



^The proportionality constant is smooth and order one. This is firstly because X S (—1,0) ev- 



erywhere, and secondly, at the point Y diverges, i.e. at ^/j, the boundary condition (A.4) guarantees 
that the proportionality constant is still order 1. 
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Figure 20: Validity of the adiabatic approximation. Solid(blue) curve is the difference 
between the true numerical result and the adiabatic approximation (normalized by the 
true value) {C/s{true) — (/s{adb)) /{C/s)(true), and the dashed(red) curve is the criteria 
I for the validity of approximation. 



for the validity of the adiabatic approximation. The regions where both functions 
become large (the region around Ac) coincide, as expected from our argumentation 
above. We also see that the approximation becomes better near the UV and the IR 
regions. 

In passing let us also note the physics features that cannot be captured by the 
adiabatic approximation. The same condition, namely that V'/V varies slowly, also 
leads to Xq = — | for the zero T theory. This means X — Xq vanishes in this 
regime. According to [|^, the gluon condensate is set to zero within this approxi- 
mation. Therefore, one cannot observe the phase transition at Tc in the adiabatic 
approximation. By the same reasoning, we learn that the adiabatic approximation 
becomes worst when X differs from Xq most, i.e. when the gluon condensate is 



largest, in other words in the region near Tc, see also figure 20 



Most equations simplify greatly with ( |C.1| ). In particular, the coefficient d{(j)) 
of ( |A.15|) vanishes in the fluctuation equation ([A.13| ) for tu = 0. Therefore with the 
same arguments of Appendix |B| we are lead to the conclusion hn = 1 in the adiabatic 
regime. Then, equation ( |3.23[ ) follows immediately. 



D. Equivalence of the axial and the 



gauges 



According to the standard AdS/CFT dictionary the metric fluctuation hu = /i22 = 
^"^^ in the axial gauge h^m = 0, whereas our computation of the bulk 
= gauge, following . 



/i33 is dual to 
viscosity is carried out 
result is independent of the gauge choice, by performing a gauge transformation 



In |53[, it is shown that the 
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between the two gauges and showing that this does not affect the couphng of the 
fluctuation to the corresponding operator to leading order near the boundary. In our 
backgrounds, this issue is slightly more subtle, due to the logarithmic corrections to 
the asymptotically AdS geometry. 



Here we shall follow the steps in and prove that indeed the gauge choice does 
not affect the coupling also in our backgrounds. The metric in the 6(f) = gauge is 
given in (|A.8|) . Asymptotic forms of the metric functions near the boundary read (in 
the A = exp(0) coordinate) , 

^(A) = ^ + Mog(M) + 0(A), (D.l) 

5(A) = -log(M) + 0(A), (D.2) 
/(A) = 1 + C(e-^/'°^(6oA)"''). (D.3) 

We want to perform a gauge transformation from the gauge I to gauge II where, 

/ : 50 = 0, Sg,, = dtag[-fe''^hoo, e^^hu, 6^^/122, e^^hs, e^^^], (D.4) 
II : 6<f) = C V 0, 6g~^, = dtag[-fe^^hoo, e^^/in, 6^^/122, e^^hs, 0]. (D.5) 

Under an infinitesimal gauge transformation (^^, the metric functions and the dilaton 
transform as 

5(j) = S(j) + C^d^ct), 5~g^, = 5g^, + V^C. + V.C/.- (D.6) 
The symmetries of the problem dictate that = = = and = (^°(t,0). 

When applied to ( p.4|) in order to get ( p.5|) these transformations reduce to the 
following equations. The dilaton equation yield 50 = and the 55, 50, 00 and 11 
components of the second equation in ( p.6|) respectively produce ||53| : 



C^' + {B' - + ^ = 0, (D.7) 

C° - -J^^' = 0' (D-8) 
hi - hu - 2AX^ = 0, (D.9) 
hoo - hoo - {2A' + - 2C°' = 0, (D.IO) 

where prime and dot denotes derivation w.r.t. and time, respectively. 

We shall assume an oscillatory form for the t-dependence of the fluctuations, 
5x{t, 0) = e~™*x(0) and use the same symbol to denote the 0-dependent piece, with 
a slight abuse of notation. One can easily see that the following computation goes 
through with no change for a more general t-dependence. Using the asymptotic 
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forms of the metric functions above, one finds an approximate solution to (p.7 ) as, 



C'(0) ^ A 



Cl 







(D.ii; 



In order to determine the asymptotic behavior, one needs to determine /155 near the 
boundary. This can be done by using the fluctuation equation for /155 ( see 



I. 

2/ 



hi 



(D.12) 



where hu is the solution to ( |3.6| ). Again, using the asymptotic forms of the met- 
ric functions above, one finds that the solution to ( |3.6| ) near boundary (with the 
boundary condition hu — t- 1) reads. 



hu ^ 1 + c'X 



w = 0, 



(D.13) 
(D.14) 



where c, d are some integration constants. Using these in ( D.12| ) one finds that 
/i55 = C(A"^^e"^) for a; 7^ and 0{\-^^~^e">^) for w = 0. Finally, substituting 
this in ( |D.11| ) gives (for w 7^ 0) , 



Cl + 0{X 



-2b- 1 



e ''0^ 



(D.15) 



We see that the leading term goes as — 1/logr whereas the sub-leading term is 
suppressed as (9(r^) as r — )■ at the boundary. Thus we can safely ignore the inho- 
mogeneous contribution in (p.ll|) and take ~ CiA. Using this and the asymptotics 
of the metric functions above in (p.8|) now gives. 



'iwt 



(D.16) 



Finally, using all the above, one solves ( p.lO|) and ( p.9|) as (stripping off the t- 
dependence). 



bo 
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e~¥^), hu = hu-'^ + 0{X-^''~'e-^). (D.17) 



The operator that is dual to = 50 is C = tr FV(4A) ([|§) . Thus, from ( PXT3| ) 
and ( p.l7[ ) we find that the fluctuation of the Lagrangian that is proportional to Ci 
is I I 

(D.18) 



A £ = -T^ + -fr 



Expanding the dilatation Ward identity, = ^^tr to leading order in A we see 
that, ( p.l8| ) vanishes. This proves that the fluctuation of the metric function hu 
couples to both in the axial gauge and in the 6(f) = gauge. 
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E. UV subtleties 



In a theory where the dilaton is non-trivial, there is a relevant question to be asked. 
In which frame is the metric asymptotically AdS? This question is void in the dual of 
A/" = 4 SYM where the dilaton is trivial but not in Improved Holographic QCD where 
the gauge coupling is a function of the holographic coordinate. In |^ by analyzing 
the structure of string higher-curvature corrections, it was shown that it is consistent 
with the equations of motion that the string-frame metric is asymptotically AdS. This 
is required in order for the background solution to have the correct structure and 
QCD perturbation theory to emerge. 

By approximating the string theory dual to QCD by a two derivative theory and 
a dilaton potential as it was proposed in ^ this property cannot be maintained. 



It is not easy to see that the only option that can be implemented in the UV is 
an asymptotically AdS metric in the Einstein frame instead. This has as a result a 
few "stray logs" in several quantities that are calculated from the world-sheet action 
(instead of the bulk effective action). One of them is the short distance inter-quark 

4 

potential calculated in |73| which is^'^ V{r) ~ (iHsIlMl., 



Other similar cases appear in this paper, in the two observables that involve 
the string world-sheet action. The first is the drag-force calculation. The effect of 
these logs appears both as the energy of the string end-point becomes asymptotically 
large {v 1) or when the temperature becomes large T — )• oo. One example is the 
ultra-relativistic diffusion time ([4.32|) that we reproduced here 



4A'2 /5 p2 \l 



4 

The logarithmic fact is due to a factor of A^a . Another example is the large T asymp- 
totics of the non-relativistic diffusion time in ( [4.33| ). It gives exactly the conformal 

4 

result modulo again a factor of A~3. 

A similar effect appears in the jet-quenching calculation in section |^. Indeed, 
in ( ^.22| ), the combination e^'^'(l — /) vanishes logarithmically in the UV instead 
of asymptoting to a constant value. This occurrence perturbs the structure of the 
Wilson loop configuration near the boundary, but as shown there does not affect the 
calculation of the jet-quenching parameter. 



A perturbative calculation at NLO of the diffusion time gives, |74| 

^ ^^'^ al [- log g + 0.07428 + 1.8869c/] (E.2) 



TpQCD 3M 



^'"^Interestingly, it was argued in Q that this fits better Quarkonium data than the CorneU 
potentiaL 
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while for large Nc N' = 4: SYM it is obtained 

1 A^T^ ' 



log ^ + 0.4304 + omiVx 

V A 



(E.3) 



TpsYM 127rM 

Such perturbative asymptotics are not visible in the NG action. 

There is however an important issue here: when and where we can trust the 
standard Nambu-Goto world-sheet action. The structure of the vacuum solution 
near the boundary, advocated in [0, suggests that since curvatures are high in 
that region, care is needed when using the NG action in that regime. On the other 
hand, the UV behavior remains qualitatively correct although in its details it may 
be revisable. 



F. The UV asymptotics of the integral (5.21) 

We now turn to estimating the integral 

4 

dr r A~3 dr 



used in section ^ when e — > 0. 
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Near r = in the Einstein frame, [3S 
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changing variables to m = — log(Ar) we obtain 

J(e) = ^ ° / du wte-" 
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We now use 



du use " = r 
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to finally obtain 
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valid as e — )■ 0. 
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